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ABSTRACT 

/ 733 ? 

There is now considerable evidence for the existence of 
an interaction between the solar wind and the geomagnetic field. 
The solar wind is a highly ionized plasma and as such may well 
influence the earth's magnetic field. It is this interaction 
that we have studied. 

We have limited ourselves to a two dimensional model to 
permit an exact solution to be obtained. The plasma wind is 
considered as impinging upon an arbitrarily oriented two dimen- 
sional magnetic dipole. We find that the field is compressed 
by the plasma and essentially confined to form a cavity in the 
plasma wind. The thickness of the boundary layer within which 
the plasma and field Interact is found to be small. We have 
determined the equation of this boundary surface and the extent 
to which the field is changed by such a confinement. 

We have also considered the stability of such a system. 

We find that there are certain regions of instability, but for 
the part of the surface facing the wind it is stable. 
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INTRODUCTION 

Recently a great deal of work has been devoted to the problem of 
the interaction of the solar wind with the earth's magnetic field. The 
solar wind consists of a totally ionized hydrogen plasma continuously 
streaming from the sun. The velocity and density in the neighborhood 
of the earth are approximately 10 cm/sec. and 100 cm J respectively. 

These values have been estimated by Biermann and others by observing the 
deflection of comet tails in the solar wind. Further evidence of its 
existence is the geomagnetic storms accompanying increased solar activity. 
The mechanism for this disturbance will be made clear in the work that 
follows . 

2 0 3 

The energy density of this plasma is l/2 Nmv = 0.83 X 10 ergs/cm 
where N is the number density of the plasma particles and m the proton 
mass. If we approximate the earth's field by a simple dipole the magnetic 
field in the equatorial plane is B = 0.35(R o /R)^ where R q is the radius of 
the earth and 0-35 the earth's field in gauss at the equator. The energy 
density of the magnetic field then is B^/8jt = (.35 R Q ^/R^)^/8n. We see 
that the energy density of the field and the energy density of the plasma 
wind are equal at approximately 9-1 earth's radii. If the plasma wind 
is to have any effect on the earth's field it should be felt at 9.1 
earth's radii and beyond. 

There is direct experimental evidence for such an interaction. Rocket 
studies 1 indicate that the earth's magnetic field behaves roughly like a 
dipole out to 12-13 earth's radii where it sharply decreases to the field 

1 A Radial Rocket Survey of the Distant Geomagnetic Field, C.P. Sonett, 
D.L. Judge, A.R. Sims, and J.M. Kelso, Space Technology Laboratories Report 
7320,2-13- 
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of Interplanetary space (5 X 10 y gauss). 

We wish to describe the nature of the Interaction between the plas- 
ma wind and the earth T s field. In particular we shall show that the wind 
compresses the field so that a cavity Is formed in the wind. This com- 
pression may be seen physically from either macroscopic or microscopic 
models . 

First let us consider the macroscopic model. The solar plasma has 
a high conductivity. Parker^* has estimated the mean free path to be of 
the order of 10^ km. Now a good conductor behaves dlamagnetically. 

Any magnetic field which exists initially In the conductor will tend to 
persist. If the field is zero Initially it will remain zero. This pheno 
menon of diamagnetism may be explained by Lenz's law. If the magnetic 
field does change there will be an Induced electric field. The Induced 
currents will flow in such a direction so as to oppose the change in the 
field. The greater the conductivity the greater the opposition to change 
Thus if we Imagine a plasma streaming toward the dipole from Infinity 
where the magnetic field is zero the magnetic field must remain zero for 
all time in the limit of infinite conductivity, i.e. the magnetic field 
in the plasma may not change. 

We may also envisage this exclusion of the field from the plasma 
from a microscopic picture. Since the mean free paths are so enormous 
a particle description is certainly better suited to describe the inter- 
action* Charged particles moving in a magnetic field themselves create 
magnetic fields. Imagine a uniform magnetic field in the negative half 


E-N. Parker, Interaction of the Solar Wind with the Geomagnetic Field, 
Phys. of Fluids, 1 , 171 (1958). 


7 




3. 


space y < 0 and zero in the positive half space y > 0. Consider a 
particle incident on the field-vacuum interface along the y-axis from 
y = + oo . it will he curved by the magnetic field moving in a semicir- 
cular path with center of curvature on the positive or negative x-axis 
depending on the sign of the charge. In either case the particle will 
decrease the field within Its orbit (i.e. within the semicircle) and 
Increase the field outside. We may now generalize to consider the prob- 
lem of a plasma Incident on the interface. For sufficiently large num- 
bers of particles the magnetic field may be materially reduced within 
the region where the particles are turning, and Increased beyond this 
region- Thus we may speak of a compression of the magnetic field by 
the plasma stream. We shall study this problem In more detail later. 

In our particle picture we may describe the Interaction as follows. 
The plasma particles move in straight lines (because of their long mean 
free paths), penetrate a short distance into the field, and are then de- 
flected out. 

A more rigorous justification of these statements will follow in 
the next section. 


DEPTH OF PENETRATION 

The problem to be described now has been considered by Dungey 1 and 
2 

Rosenbluth . Our reason for Including It here is partly for complete- 
ness and partly because we feel that an important feature of the problem 


J.W. Dungey, Cosmic Electrodynamics (Cambridge University Press, New York 
1953), Sec. 8.3. 

M. Rosenbluth, In Magnetohydrodynamics, edited by R. Landshoff (Stanford 
University Press, Stanford, California, 1957)* 


8 




4 . 


has been Ignored. 

We would like to know how far the plasma will penetrate Into the 
magnetic field. To this end we shall Investigate the following two 
dimensional problem. The magnetic field Is constant in the half space 
x < 0 and approaches 0 as x approaches Infinity. The field Is every- 
where parallel to the z-axis. The plasma is Incident on the field from 
x =s + ». We assume that all the particles are moving parallel to the 
x-axIs before encountering the field. We expect the path of the par- 
ticles to be curved as they enter the field and be deflected out. The 
protons having the greater momentum will penetrate further than the elec- 
trons, thereby building up a space charge. The electric field of this 
space charge will tend to pull the electrons in with the protons. The 
paths of the particles are roughly described in Fig. 1. We observe that 
the particle current is responsible for the decay of the magnetic field 
in the x-direction. The discontinuity in the magnetic field between 
x = 0 and x = » is equal to the total particle current in the y-direction. 

The solution of this problem is determined by a self-consistent 
solution of the equations of motions of the particles and Maxwell's equa- 


tions. The boundary conditions are that B = B q at x = 0 and B 0 as 
x oo. The electric field must be 0 at x = 0 and approach zero as 
x -► + oo. Let (u^x), v (x), 0) be the (x,y,z) components of the outgo- 
ing protons. From symmetry we see that the velocity components for the 


Incoming protons are (-u^,v ,0) . The corresponding velocity components 
of the electrons are (u^v^, 0) and (-u^, v e ^)“ Let n p( x ) n e ( x ) 
the number density of the incoming and outgoing protons and electrons 


respectively. The magnetic field is represented by a vector potential 
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having only a y-component A(x), and the electric field hy a scalar poten- 
tial ^(x). The number density and velocity at x = °o are N and U. 

We may integrate the y-component of the equations of motion 


v 

P 



v 

e 


eA 

m c 
e 


0 , 


( 1 ) 


The energy equation gives 


2 

u + v 
P P 


2 2ed 2 2 2ed TT 2 

+ — ^ = u + v = U 


m 


m 


( 2 ) 


Particle conservation requires 


n u =s n u = NU 
p p e e 


Maxwell’s equations are 


d A -8 rce f N 

cbc 2 c PP e e 


and 


^-4 = -8jte(n - n ) . 


dx 


P e 


( 3 ) 

(M 

( 5 ) 


(The currents and charge densities have been doubled to include both in- 
coming and outgoing particles.) 

The above system of equations apply in region II only. If the elec- 
tron terms are dropped we obtain the equations appropriate to region I. 

It goes without saying that these equations cannot be solved exactly. 
We can however find a valuable first integral, (it is at this point that 
we depart from Dungey’s analysis.) 
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These six simultaneous equations can be reduced to two simultaneous 
equations relating the vector and scalar potentials. They are 



We may simplify the right hand sides 



If we multiply the first equation by -dA/dx and the second by d^/dx 
and add them, the right-hand side becomes an exact differential: 




7 . 


Integrating ve obtain 


2 2 

E H 

Bit Bit 


2NU(m u + m u ) 4* constant. 
' p p e e y 


To evaluate the constant ve let x « where E = H = 0 and u = u = U. 

P e 

Therefore 


~ 2 H 2 


I? “ £ - 2NU ^p U p + “eV ’ ^“p + m e ) ' 


( 6 ) 


At the interface (x = 0) E = 0, H - H , and u = u =0. Hence 
N o' p e 


H 2 = l6jtNU 2 (ra + m ). 
o v p e 


Thus ve have found vhat magnetic field strength is required to hold back 
the plasma without obtaining a complete solution to the problem. 

We might note that the integral obtained In equation (6) could be 
found more Bimply using the Maxwell stress tensor. Consider the volume 
of plasma contained in a parallelpiped of unit cross sectional area. The 
face a is located a distance x from the plasma-field interface, and the 
face b at *>• (See Fig. 2.) The net force acting on the plasma in this 
volume must equal the time rate of change of the particles within. The 
force on face a Is H /On - E /8n. The force on face b is zero since the 
field vanishes at <». The forces on the remaining four faces cancel by 
pairs. We shall equate this force to the time rate of change of the mo- 
mentum of the particles. In the initial configuration the total x-com- 
ponent of momentum is zero since there are as many particles moving to 
the right as to the left. The momentum change which occurs in the time 
dt will then be equal to the total momentum of the system after a lapse 


? 2 


of time dt 
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(We have omitted the motion In the y-direction as this would only compli- 
cate the picture and is of no interest in calculating the change in the 
x-component of momentum.) Those particles which were at the surface a 
and moving to the left are at the surface after a time dt. Those 
moving to the right are located at / . V and b" are defined in a simi- 
lar way . The change in momentum then is the total momentum of those 
particles between a' and b" which were between a and b at time t . There 
will be other particles between a' and b" which moved in during the time 
dt. These must not be counted. The momentum of the particles between 
a" and b' is zero. The momentum of the particles between a' and a M is 

-2u dtn m u and the momentum of those between b' and b" , 2UdtNm U. 

P P P P p 

We must add similar expressions for the electrons. Equating F and dp/dt, 


2 ? 

E h 2 2 2 

-a 73 — = 2n m u +2nmu - 2N(m + m )U , 

ojt on PPP eee v p e' 2 * * * * * * 9 


but since n u 
P P 


n u = NU 
e e 


2 2 

75-7 - 7 T— = 2 NU(m u + m u ) 

cm cm P P e e' 


2NU (m + m ) 
P e 


which we see is eq. ( 6 ). 

It is not possible to proceed any further with an exact solution. 

We would however like to estimate the rate of decay of the magnetic field, 

i.e., how far from the vacuum-plasma Interface must we go before the mag- 
netic field becomes a small fraction of the vacuum field. It will be 
important in our subsequent work that this distance is small compared 
to the dimensions of the cavity. 

Pungey has presented an approximate solution to this problem. He 
assumes that the electron density is everywhere approximately equal to 
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the proton density. He has justified the approximation using the full 
set of equations (1,2, 5)« In <foing so he has overlooked the possi- 
bility of an absolute change separation such as would occur in region I. 
In this region the electron terms must be omitted from the equations and 
hence are not to be included in an estimate of the relative densities. 

We will first find an upper limit on the thickness of region I and 
then verify that Dungey’s solution is valid in region II. Having stud- 
ied both regions we may then estimate the thickness of the current 
sheath. 

We return to the integral we found in eq. (6). 

E 2 = H 2 + l6jrNU(m u + m u ) - l6rtNU 2 (m + m ). (7) 

v p p e e v p e v 

Now at the interface between region I and II 


u < U 
P 

H < H 

o 


and 


E = 8*e 


dx » Snte 


u 

/ 


N dx = 8neNx,_ 


o o 

where x q Is the distance between the plasma-vacuum Interface and the 
Interface between I and II. Substituting into eq. (7) we find 


(8jteNx ) 2 « E 2 = H 2 + l6*NUm u - 2NU 2 (m + m ) < 
o' p p P e 

< H 2 + l6jtNU 2 m - l6jtNU 2 (m + mj 
o p P e 

= l6*NU 2 m . 

P 
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So that 


or 


2 


x 

o 


« 


m U 2 
^jte 2 N 


( 8 ) 


x « 8 meters, 
o 

Where 8 m. is a large distance in laboratory experiments, it is quite 

7 8 

small in comparison with the dimensions of the cavity, 10 - 10 m. 

We might obtain a better physical insight into this upper limit on 
the charge separation by considering how much energy is available in the 
plasma stream for charge separation. Let us imagine a block of plasma 
of unit cross sectional area and width a moving to the left with a velo- 
city U. The plasma block now runs into an army of Maxwell demons who 
instantaneously stop all the electrons and allow the protons to pass. 

We shall simplify the calculation by assuming that the protons move as 
a block, and ve ask how far this block will move before coming to rest. 
When the protons are at rest the energy stored in the electric field 
must equal the initial kinetic energy. If the protons are displaced 
a distance x q (see fig. 3.) from the electrons the electric field is 


E = ^neNx . 

o 

Equating the field energy to the kinetic energy 

2 2 *r 2 2 

p 16* e N x p 

E l On = a = l/2 Nm p U 

so that tt 2 

2 m p U 

4jtNe 2 

which is identical to the upper limit for the change separation obtained 
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in eq. (8). (We have tacitly assumed that h » x q .) The fact that the 
numerical coefficient is the same is of course fortuitous. 

Nov let us consider region II. Here ve shall show that Dungey’s 
solution is quite good. He assumes that the electron and ion charge 
distributions are approximately equal. This assumption is to he tested 
a posteriori. From eq. ( 3 ) ve see that 


u = u . 

e p 


The y-components of the electron and proton velocities are related by 
eq. (1) 


m v = -m v 
P P e e 


so that |v | « |v j. We may therefore neglect the proton current in 
p e 

eq. (k) so that 


d 2 A 8jteNU V e 


dx 


c u 


( 9 ) 


If we now eliminate <f> from eq. (2) and substitute for v and v from 

p e 

eq. (l) ve find that 


2 2 

2 e A tt 2 

mu + tt = m U 

P P 2 p 

r * me 

e 


( 10 ) 


neglecting terms of order m /m . We now solve eq. (10) for u and eq. (l) 

e p p 

for v and substitute in eq. (9) remembering that u = u , 
e P e 


a A 8*eNU 

7 ^ * c 

dx 


eA 

m c 
e 




(id 


m m c 
e p 


1 ^ 
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a = k/k 


where 


mm c 2 !^ 

A 2 _ JU2 

o 2 

e 


X 1 = x/x 

' r 


2 m e C 

x = ? 

° 3*Ne‘ 


Eq. (11) then becomes 


f l - a 


Clearly a ^ 1 so we set a = sin 9. Integrating 

x 1 = - tn(tan ©A) - 2 cos(9/2) + tn( ~J2. - 1) + y/2 


For small 9 


in g/4 


so that 


9 = be' 


H * 30 e 


The field then will decrease by a factor l/e when x 1 =* 1, i.e. 



& 500 m. 


for N = 100 cnT^. This is again much less than the dimensions of the cavity. 
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We must now justify the assumption made earlier that the electron 
and proton densities are approximately equal. We may do this hy calcu- 
lating 

~ -8ite(n - n ). 
dx P 


Now 


, m e , 2 2 2> 

* ’ 2i (u e - V e ' U } 


and from eq. (10) 


2 it 2 
u = U - 

e 


2 2 2 

= u 2 - u 2 S 

m m c A 

e p 0 


U 2 COS 2 © 


Also 


0 2,2 m „ .2 m 0 

v 2 - - — = — E U 2 — — = — E— u 2 
v e 22 m . 2 2 U 

me e A m 
e o e 


sin © 


( 12 ) 


so that 



(13) 


e 

We observe that if m 0 were to equal m^ the potential would be 
zero and consequently no charge separation vould occur. This is to be 
expected. If the particles have the same mass they must penetrate to 
equal depths. 

Differentiating eq. (13) twice ve find 


d 2 d( , x 8*NeU 2 

72 - (m P “ V 5" 

dx y in c 


2-2 cos 9 - 


2 

sin Q 
cos 0 




so that 


n - n a 
e p 


( m p~ m e )Nu r 


m c 
e 


2-2 cos 0 - 


2 

sin 0 
cos 0 


]• 
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We see that the plasma is far from neutrality at 9 = n/2, i.e. at the 

point where u = 0 (n = NU/u -* »). Although n and n differ con- 
e e e s p 

siderably near x ■ 0 (9 = jt/2) the approximation may still be a good 

one. We must go back to our calculation of the field decay and investi- 
gate what is a reasonable measure of large deviation from neutr ali ty. 

In calculating the field from eq. (9) we substituted for u the 

e 

value obtained from eq. (10) for u assuming that u = u . Let us see 

P p e 

how far n^ and n must differ before this becomes a poor approximation. 

Let n - n =5. Now 
p e 

n u = n u 
p p e e 


or 


n 

u = — £ = 
e n 

e 


1 + S- 

s n e 


We require that 


n - n 
P e 


« 1 . 


Now 


since 


NU „/ 

n e = — ^ U/cos 9 

e 


u = U cos 9 
e 


from eq. (12). So that 


n - n (m - m )U* 

'n ' -rs o ' 


e v v e' , 2 P 

= — * g ( 2 cos e - 2 cos 0 - sin 0) 


n 


m c 
e 


m _ 

- JL_ 


u 2 


m c 
e 


(2 cos 9 - cos ^9 - 1) < .002 or .2$. 
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This completes our discussion of the rate of decay of the field in 
the stream. We have shown that the region where the field differs appre- 
ciably from zero is much smaller than the dimensions of the cavity. (We 
shall, in the future, refer to this region as the current sheath.) Since 
the mean free path is so large the plasma constituents will behave as 
free particles until they reach the cavity boundary where they are reflec- 
ted. We shall now consider the problem of the reflection of particles 
with an arbitrary angle of incidence. 


REFLECTION FROM THE CAVITY BOUNDARY 
We wish to show now that if the current sheath is small the particles 
will be specularly reflected from the sheath. We shall consider charged 
particles Incident on a magnetic field as above except their velocity com- 
ponents at infinity are now (U, V, W). The magnetic and electric fields, 

B = B(z) x and E = E(z) z, approach zero as z -*«■ We shall show that 
after reflection the velocity components will be (U, V, -W) so that the 
angle of incidence is equal to the angle of reflection. (See Fig. 4.) 

The fields may be expressed in terms of a vector and scalar potential 

A = A (z) and ^ = j^(z). The path of the particle is determined by the 
«/ 

conservation lavs. 


u * U 

v + ^(j4 = y 

me 

2 . 2 eef „2 „2 

v + w + — £ » V + W 
me 


w 

(15) 

(16) 


where (u, v, w) are the velocity components at some intermediate point of 
the motion. The constants of the motion in eqs. (l4, 15, 16) have been 
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evaluated at t = - ». We now look at t = + ». Assuming that the field 
has sufficient depth that the particle is reflected and not transmitted 
then z-^ + oo as t -► + «>. Eq. (15) requires then that v V as t -* + »• 
Substituting u = U and v = V in eq. (l 6 ) and observing that ^(00) = 0 


2 

W s 


W 


2 


so that 


w 


t=-H» 


« + w 


V 


trs-00 


- w . 


We might generalize this conclusion somewhat to include curved fields 
with the proper symmetry. However, for our purposes it will be sufficient 
to require that the field may be approximated locally by a plane field 
whose principle gradient is perpendicular to the field (i.e. within the 
orbit of a particle, the field does not vary appreciable along the field 
lines) . 


DEFINITION OF THE PROBLEM 

Let. us now return to the problem of physical interest, a plasma wind 
blowing on the magnetic field of a dipole. As we have demonstrated , the 
plasma constituents will move in straight lines up to the cavity boundary 
where they are specularly reflected. The region of interaction (current 
sheath) separates the plasma from the field. Since the sheath thickness 
is small we may treat it as a boundary surface separating the two domains. 
In the plasma domain we have streaming particles making elastic collisions 
with the boundary and in the field domain we have a static magnetic field 
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confined to a cavity by the boundary surface. The equilibrium shape of 
the cavity surface will be determined by a condition of pressure balance. 
The magnetic pressure within the cavity must balance the pressure of the 
plasma wind from without. We may then formulate our Idealized problem 
as follows : 


Vxi> 0 


(17) 


everywhere inside the cavity except at the dipole where 
B - M/r^ sin 9 0 - 2M/r^ cos © r 

as r -+ 0. 

V ♦ B = 0 

inside the cavity, 

B = 0 
n 

on the surface, and 

B^/8?t = 2NmU^ cosSc 


(18) 


(19) 


( 20 ) 


on the surface. Here X is the angle between the normal to the surface and 
the direction of the incident wind. 

In eq. (17) ve neglect all currents within the cavity such as atmos- 
pheric currents and currents in the Van Allen belts, or else we lump them 
in with the dipole field. The boundary condition (eq. (19)) is imposed 
since the field is zero outside the cavity and the normal component must 
be continuous. Eq. (20) expresses the condition of pressure balance, i.e. 
we are looking for a steady state solution. 

This problem differs in an essential way from the usual boundary 
value problem. If the magnetic field were expressed as the gradient of a 


r, 
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scalar field, the problem -would be determined by finding a scalar solution 
to Laplace’s equation having the appropriate singularity at the origin and 
satisfying two boundary conditions . Normally one is given the boundary and 
one boundary condition. Here the second boundary condition must determine 
the boundary . 

Unfortunately, even with these simplifications, the problem is still 
much too difficult to handle. The existing theoretical work on the shape 
of the cavity boundary in three dimensions consists entirely of free hand 
drawings. There is one exception. David Beard^ has presented an approxi- 
mate solution to the problem. The validity of the approximation is ques- 
tionable and there exists no test of its accuracy. We shall present here 
the exact solution of two analogous two dimensional problems, namely, the 
interaction of a plasma wind with the magnetic field of an infinite line 
current and with the field of two lines carrying currents in opposite 
directions (i.e. a two dimensional dipole). We do not propose that these 
two dimensional problems will approximate the three dimensional problem. 

We present these problems firstly because they can be solved exactly; 
secondly, they are of interest in themselves, giving us physical insight 
into the nature of the interaction between plasmas and fields; thirdly 
they can be approximated experimentally so that the theory may be tested; 
and lastly they will serve as a test for approximate theoretical solutions 
to the three dimensional problem. Regarding this last category, we shall 
apply Beard’s approximation method to the two dimensional problems and 
demonstrate explicitly the validity of the method. 

1 David B. Beard, The Interaction of the Terrestrial Magnetic Field with 
the Solar Corpuscular Radiation, J. Geophys. Research, 65 , 3559* I960. 
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INTERACTION BETWEEN A PLASMA WIND AND THE 
FIELD OF AN INFINITE LINE CURRENT 
We have developed two methods for the solution of the problem-type 
defined by eqs. (17,18,19, 20) . The first, while offering considerable 
advantages In algebraic simplicity, is restricted in the class of prob- 
lems It can handle. In particular it may not be used on the two dimen- 
sional dipole. To present this first method we consider a plasma wind 
Incident on the field of an infinite line current. (See Fig. 5 .) The 
plasma velocity is everywhere perpendicular to the current. The field 
Inside the cavity boundary must satisfy Maxwell's equations 
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We now reformulate the problem in complex notation. Since 
2 SgCftal , 7 • B 4 i |V X B| . 0, 

B must be a function of z* where 

B =» B + IB . 

x y 

The singularity at the wire required that 
B 2Ii/z* as z -► 0. 

To formulate the boundary conditions consider 
B*dz = B • ds + i | B X ds | . 

The imaginary part must be zero (tangency condition). The real part is 

i 

( 8 rtP 0 )^ cos © ds 

from the condition of pressure balance. Since cos G » dy/ds we may 

state both boundary conditions in the single equation 

l 

B*dz = ( 87 tp o ) 2 dy. 

Our problem then is to find a complex function B such that 


B = B(z*) 

( 21 ) 

B -*• 2Ii/z* os z *♦ 0 

( 22 ) 

B*dz = ( 8 xp o )^ dy 

(23) 


on the surface. 

We now use the unknown field B to generate a conformal mapping. 

Let 

w(z) - l/iB* 

1 J.D. Cole and J.H. Ruth, Physics of Fluids, 2, 62^ (1959)* 
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where 

w = u + Iv. 

(The choice of the mapping function is determined by the simplicity of 
the resulting boundary value problem.) Putting this in eq. (2j) 


dz = (8rtp Q ) 2 dy iw 


(24) 


on the surface. Taking the imaginary part 

U = (8kp o ) -? = U Q , 

a constant. We see that in the w plane the problem has been reduced to 
one of known boundary and is thereby considerably simplified. 

The other boundary condition contained in eq. (24) might be obtained 
by taking the real part. Instead we proceed by considering the quantity 

- p- = i(8jtp )* §£ 

v dw o dv 

on the surface. But on the surface dv = i dv and Is pure imaginary. There- 
fore 

Im(l/v dz/dv) = 0 


on the surface.- i.e. when w = u 4* Iv. 

o 

The singularity condition near the current is that 
B -*■ 2Ii/z* . 

Therefore 

z 2Iw 


and 

l/w dz/dw -► 2l/w as w -► 0. 
These condisions will be satisfied by taking 
l/w dz/dw = 2i(l/w - l/w Q ) 
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vhere (see Fig. 6) 


v = v - 2u . 
o o 

We find then for z 

z =* -4 1 u In 
o 

The constant of integration a is determined by the singularity condition. 
Let a = ~2u . Then 



z » -^Iu q ln(l - v/2u o ) (25) 

and approaches 21 v as it should. 

If ve evaluate v in eq. (25) on the surface (v = u q + iv) ve deter- 
mine the surface 


so that 


/ u + iv^ 

= -4lu In (l 

o V 2m 


, 2 2 \i 

(u + v ) 2 

x = -4lu in 

s o 


2n 


( 26 ) 


(27) 


4lu tan v/u . 
o ' o 


(28) 


vhich are parametric equations for the boundary surface. These results 
eire plotted in Fig. 7- 

Eq. (25), besides giving us the shape of the cavity boundary, can 

be inverted to give us also the magnetic field. Solving for w 

Z' -z/4lu N 
v=2u o M- e °) • 


But 


v 


1 

IB* * 
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Therefore 




-xAlu 

e sin yAlu 

o 


x 2 - 2 xAlu 2 

cos y/4lu Q ) + e sin yAlu Q 


and 


B 


1 

2u 

o 



1 - cos yAlu 
' o 


-x/4lu X 2 -2x/4lu Q ^ 

e cos y/4lu J + e sin y/4lu Q 


We observe th at d oes not approach zero as y approaches zero hut 
Instead lAu or v . Thus the field Is Increased on the windward 

side of the line current and decreased on the leeward side. One might 
explain this phenomenon by Imagining that the cavity boundary has com- 
pressed the field lines between the cavity and the line thus strengthen- 
ing the field. The field on the windward side points generally In the 
positive y-directlon. We expect then that the surface currents on the 
boundary will produce a field In the y-direction near the line current. 


INTERACTION BETWEEN A PLASMA WIND AND THE MAGNETIC 
FIELD OF A TWO DIMENSIONAL DIPOLE 

This problem is formulated In exactly the same way as the previous 
problem except that now the singularity Is that of a dipole Instead of 
a line current. We take the dipole to be perpendicular to the plasma 
flow. This change alters the problem radically. It Is not possible to 
use the previous technique here.. The reason will be made clear shortly. 
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We begin again with the relevant equations, 
V x B - 0 

everywhere except at the dipole where 



as r -► 0. Here M is the dipole moment. 

V • B » 0 
everywhere . 

The field must be tangent to the cavity surface so that 
B 

JL = & 

B dx 

x 

and the magnetic pressure must balance the pressure of the stream, 

E (\ 2 + B /) “ Po C ° s2 * 

where 

p « 2Nmv 2 . 
o 

Again we find a considerable simplification by expressing these 
conditions in complex notation. The two Maxwell's equations become 

B = B(z*) (29) 

and the two boundary conditions 

B* dz = B • ds + i B X ds = (+) ^/8jtp o dy. (30) 

The singularity at the dipole is 

B* -+ — i as z 0. (31) 

z 
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The reason for the (+) in eq. (50) can he seen from the geometry of 
the field. We have made a rough sketch of this in Figure 8. 

If ds is a vector tangent to the surface and the sense is defined 
in the usual way (moving along the surface we keep the internal domain 
on our left) ve see that B*ds (i.e. the real part of B* dz) is negative 
between the two neutral points (n) and positive outside. 

The magnetic field reverses direction at the neutral points and is 
therefore zero there. Since the magnetic pressure Is zero the particle 
pressure, and hence the slope, must also be zero. 

It Is this reversal of sign that prevents us from using the previous 
method of attack on this problem. If ve again used the field to define 
a transformation, the transformed boundary shape would no longer be a 
plane. This plane was defined by the condition that u = l/ = u q . 

This condition becomes in the present case u = + u q depending on whether 
we are inside or outside of the neutral points and the transformed geometry 
would become 


u 


Fig. 9* 
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The geometry is complicated further hy the fact that the boundary 

lies on a double sheeted Riemanian surface, (This follovs from the fact 
2 

that v -► z as z 0, so that where an argument of 2jr described the z- 
plane an argument of kjc is needed to describe the v-plane.) 

For these reasons we must develop a new method of solution. In 
brief we may describe this method as follows. We assume there exists a 
mapping w = w(z) which transforms the boundary onto a circle. The boun- 
dary conditions then define a set of requirements which this transforma- 
tion must satisfy. We may solve these equations for v(z). Then knowing 
the transformation we can map the circle In the v-plane back onto the 
z-plane to find the boundary shape in the physical space. 


AMLYSIS 

We introduce a potential function $ such that 
d<J> 


B* =* 


dz * 


(32) 


The conditions required of $ are that 
$ * *(z). 

This follows from eq. (29). The singularity at the origin requires that 


. 2M 

S> -► — as z -► 0. 
Iz 


And finally the boundary conditions become 


(33) 


H dz = (+) d y 


W 


We now assume that there exists a complex transformation w = w(z) 
which maps the unknown boundary In the z-plane onto a circle of radius 
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a in the w-plane. We further require of the mapping that it be an iden- 
tity transformation near the origin, i.e. w -► z as z -* 0. We shall see 
that this added restriction does not overdetermine the problem. There 
are many transformations which map the surface onto a circle and of these 

we choose that one which will preserve as nearly as possible the original 

2 

geometry. If for example w -► z near the origin we would be dealing with 
a double sheeted surface. 

Our problem has now reduced to finding a function <J>(z(w)) = 0(w) 
such that 

$ -► ~ as w -► 0 (35) 

Iw 

and 

d$ = (+) (36) 

on the surface of a circle of radius a. It must be recalled that the (-) 
in equation (36) is to be applied at points on the boundary within the 
neutral points and the (+) to points outside the neutral points. Now 
the transformation will determine where these points are mapped on the 
w-plane. We shall place one further restriction on the transformation. 

We require that the neutral points are mapped to the top and bottom of 
the circle (i.e. ^ = x/2 and </> = -tc/ 2) . We have now placed a number of 
restrictions on the transformation. It may be that there are too many 
and the transformation is overdetermined. We shall show that It is not 
by exhibiting a transformation which satisfies all these requirements. 

From eq. ( 56 ) we see that Im(d$) = 0 on the boundary, i.e. on the 
circle. This Implies that 

Im $ = constant on the circle. (37) 
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We can see this condition physically in this- way. Let $ = <{> + ity. 


Then 


B* 


.^ + i4 

” ’ __ 1 -T _ 


or 


dz dz 

Since $( z) is analytic ve may set dz = dx so that 


B* = + i ^ 

OX 


B = 

X ox 


and 


--s 


chlr ck 


where ve have used the Cauchy-Riemann equations. Therefore 

B = Be+Be=xV0. 

xx y y 

The magnetic field lines then are the family of orthogonal trajectories to 
the equipotential surfaces ^ = constant. But we know that since f and \jy 
are conjugate functions that ijr = constant are the family of lines orthogo- 
nal to ^ = constant. Therefore the magnetic field lines are coincident 
with the family of lines ^ = constant. Now on the cavity surface the mag- 
netic field lines must he tangential so that the surface is itself a 
field line. Therefore \jr « Im $ « constant on the surface. Ve shall choose 
this constant to he zero. 

We shall attempt now to construct a trial solution for the potential. 


Let 


2M 

iw 


2M 

ia c 


v 


( 58 ) 


This choice clearly fulfills the requirement of eq. ( 35 ) since w -► z as 


0. It also satisfies eq. ( 37 ) since on the circle w = a e 




and 
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. 2M 2M i.i 4m . / 

ia ia a ^ 

so that Im $ = 0 on the circle* 

There remains only to investigate the real part of eq. (%), i.e. 

- — cos^d^ = + «yS5tp dy. 
a o 

This equation is readily integrated, 


y' = 


- 2 - sin <ji 
< sin 

2 - sin (f> 


- it < <f> < - jt/2 

- jt/2 < < jt/2 

tc/ 2 < ^ < jt 


where 

a - v /Sit p Q 

y ' = 5 m y ' 


We shall find it more convenient to write this equation in this form 



( 50 ) 


We wish now to find a function z' = z'(w) such that z' is propor- 
tional to w near the origin and whose imaginary part satisfies eq. (4o) 
on the boundary. This function will define the transformation from the 
z-plane to the w-plane . 
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We shall construct this function piecemeal. Let us "break up the 
boundary value problem into tvo parts 


2 - sin <j> 


-2 - sin 


sin <j> = 2 


0> + sin ^ <-l 


Consider the function 


1 . 1 
F 1 = 2 ^ ln 2 2 • 
w 2 


It Is not difficult to show that 


Im F 1 ^ 


when evaluated on the boundary. (The variables are defined in Fig. 10.) 


Also, If 


1 W 2 

F„ = - In -% +21 

d * v d 

3 


Im F 2 = -1 


when evaluated on the boundary. 
Finally, let 

Fj = i/2(a/v - w/a) 


Fj = sin ^ 


on the circle. 
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Let us now evaluate the imaginary part of 

F,F 0 + 2F.. + real constant 

J2 1 

on the boundary* We see that it is Just y 1 . (The addition of a real 
constant will not affect this result*) We therefore set 



We have chosen the constant so that z' -*• 0 as w 0* In particular 



jra 


as w -► 0, or 


16M 

7ta P 


v 


as 


Nov we assumed that z ■» v as v -* 0. 
so that the coefficient of v is one, 


v -► 0 • 

We therefore choose the radius a 
i .e. 


2 16m 

a 38 ~ ' ' 

n ^/8jrp o 

The real part of eq. (4l) evaluated on the boundary is 


x' = - tn jQL±_£2S_iif 
cos 



+ sin tn 



(42) 


( 43 ) 


Eqs. (39) and (43) are parametric equation for the cavity boundary. 
This has been plotted in Fig. 12. 
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There are several interesting features of this shape. First we 
notice that the neutral points lie to the windward side of the dipole. 

We believe that this result may be carried over to the three dimensional 
problem. We have considered several simplified cavity shapes and in all 
cases the neutral points lay on the side of greatest confinement. As 
a simple example consider a dipole field confined by two planes intersecting 
at right angles as shown in Fig. 13. We wish to show that the neutral 
points lie to the right of the dipole. The boundary conditions are that 
the magnetic field must be everywhere tangential to the surface. This 
condition is satisfied by placing three image dipoles as indicated in the 
figure. If these are two dimensional dipoles the field lines are circles. 
There is one field line which is common to all four dipoles and this is 
the circle with center at the right angle. It is clear that the field 
is zero at the two points where the circle intersects the planes. These 
are the neutral points and again they lie on that side of the dipole 
where the field confined by the planes the neutral points would be altered 
slightly but would still lie to the right of the dipole. 

Another interesting feature is the distance from the nose of the 
cavity to the dipole, 1.35 VW • This distance is often est ima ted 

by locating that point on the earth- sun line where the energy density of 
the magnetic field of the dipole in free space would equal the energy 
density of the plasma wind. Where the two energies are comparable we 
expect that the plasma will interact strongly with the field. We equate 
then 
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n 2 'Q 1 „ 2 

B / ojc = ^ Nmv 

with 



r 


and find for the cut off distance 


r « 1.41 % '2M/ V8«P 0 . 

We see that the estimate is good to within 

The magnetic field within the cavity may be determined from 
eqs. (32,38,4l). It is of particular interest to investigate the field 
near the dipole. We find that at the commencement of a magnetic storm 
(i.e. at a time of increased solar activity) the magnetic field at the 
surface of the earth at latitudes below the auroral zone is increased. 
We shall now show that the magnetic field due to the currents in the 
current sheath of the cavity boundary will enhance the field of the 
dipole in the equitorial zone. 

If we expand eq. (4l) in a power series we find 


1 wf 1 v 3 
V " 3 a + 3 + 


Inverting this series we obtain 


w 


z 
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or 



We see that we have a constant field superimposed on the dipole near the 
origin. This field is 

By = - 2M 1 ^ r . 

9a^ 

such a field will enhance the field in the equitorial region and diminish 
the field at the poles (see Fig. l4). 

We should point out here that we have neglected the temperature of 
the plasma wind. Certainly this assumption is valid when the slope of 
the cavity wall is far from zero. The pressure of the wind on the boundary 
is determined by the normal component of the plasma velocity. When the 
slope is small the normal component of the velocity of mass flow may be 
comparable with the thermal velocities. We expect then that the leeward 
side of the cavity boundary should be pinched down when the thermal pressure 
becomes comparable with the magnetic pressure. 

ARBITRARY ANGLE OF INCIDENCE 

Certainly the earth's dipole is not in general perpendicular to 
the solar wind. It is of Interest therefore to investigate the effect 
of orientation on the cavity surface. 
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We might also mention at this point that the problem discussed 
in the preceding section has been solved independently by Dungey^. 

Our solutions agree. The methods are however different, Dungey points 
out that his method cannot be used except for the particular symmetry of 
perpendicular incidence. There is one other orientation which possesses 
a similar symmetry; the dipole pointing directly into the wind. But even 
here his method is not adequate to handle the problem. We shall show now 
that by using the technique developed in the previous section that we may 
solve the problem for arbitrary orientation. 

We assume that the dipole makes an angle a with the vertical. The 
procedure is exactly as before. The magnetic field near the origin 

-ia 

_ 2Mie 
B* ^ 

z 

so that 


The boundary condition on $ is that 
d$ = + «/8 itp Q dy 

on the surface. Again we assume that there exists a transformation 
w » w(z) that maps the cavity surface onto a circle of radius a. We 
also require that the domain close to the dipole is not distorted, i.e. 

w z as z -+ 0 


J. W. Bungey, Penn. State Univ., Scientific Report Wo. 135, duly i960. 
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We take as a trial solution 



This satisfies the required singularity condition near the dipole and 
the imaginary part of the boundary condition since 

$ _ sin($z!) + a) 

on the surface. It is clear that such a choice for the potential requires 
that the neutral points map onto the points ^ ^ - a and (f = - - a, 

since at these points ^ is zero, so that the field is zero, (see Fig. 15) 
Nov the real part of the boundary condition is 


y = ± 


4M 




sin(^I + a) + const. 


on the surface. We choose to express this equation in this way 


-sin(^+a) + B 


y l = 


-sin(^-fcx) + C 



sin(^+a) + A 


where we have set 
a 

y ' =— • 

Now since y f is continuous at N and N T 

B = A + 2 and C « A - 2 
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so that 


y’ = 


-sin(?4a) + A + 2 


sin(^+a) + A - 2 



sin(^+a) + A 


We cannot determine A except to say that A = 0 when a 0 and A = 1 
when a * rt/2. We shall leave A arbitrary and determine it later from 
the condition that z -** w as w -+ 0. 

Also it is not clear where the point x '= -« will map on the circle. 
It is obvious that when a * 0 or jc/ 2, x = -a> will map onto <f> = ir. We 
shall assume that this is true of all dipole orientations. We shall 
show that this assumption is compatible with the requirement of the 
transformation. We might instead leave it arbitrary, letting x = 
map onto <f> - jr - (3. We suspect that the condition at the origin would 
become 

iP 

z -► we 

so that (3 ss 0. However this would complicate the analysis considerably 
and in as much as a solution exists for p = 0 we may infer from the 
uniqueness of the solution that p must be zero. 

We then solve this boundary value problem for y as before and find 


that 



2ia 

n 




V 1 V 2 


+ iA + c . 


( 44 ) 
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where w^ w 2 and w^ are defined in Fig. 1 6 . 
Expanding about the origin we find 


, iE . £ + £i2 + ja + 


*a * ir 


so that A = — — and c = ^ . Then 

jf 2 


4m 


V 8rt P r 


i6m 


rta^ V8«P 


w 


so that 


2 16M 

a = 


jt V8np 


as before* 

Evaluating eq. (44) on the circle we find the parametric equations 
for the cavity boundary. 


I *: - -f 1 - *» (fS$?b) - - 1 


and 


-sin(^+a) + 2 


2a 


y s 1=1 -sin(^4a) - 2 - — 



sin(^+a) - 


2a 
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When a a 0 these equations reduce to eqs. (39) and (^3) as they must. 
We have plotted these equations in Figs, 17 and 18 for a = 37° and 
a « 90°. 


AN APPROXIMATION METHOD 

Recently Beard 1 has presented an approximate solution to the 
problem of the interaction between the scalar wind and the earth’s 
dipole field. As a test of the method we shall apply it to our two 
dimensional problem and compare it with the above solution. 

The essential feature of Beard's method is the way in which he 
approximates the field. Imagine a plane surface brought in from 
infinity compressing the field in front of it. (see Fig. 19) This 
problem Is easily solved by placing an image dipole at 0 T . If the 
plane Is Inclined at an angle a with the dipole at 0 then the Image 
dipole at O’ is inclined at an angle 2a. It is clear that the resultant 
magnetic field at the plane is just twice the tangential component of 
a single dipole. 

Beard assumes that this will be approximately true for a curved 
surface as well. The dipole field is 

2M sin 0 a 2M cos 0 a 

B = 2~~ e r 2 e 6 

r r 


1 


Op. cit. 
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(The geometry Is defined In Fig. 20.) The tangential component of the 
field then is 

2M fln \ 

2 cos (t + 0) 

r 


The pressure balance condition at the boundary is 


B t 

= P o cos X 


where 


— Lm 

B t = + — c os(\|r + 0) 
r 


so that 


cos X = cos(e - t) = + ^ c - os (* + e ) 

-y/Sxrp^^r 2 


(^ 5 ) 


The upper sign is to be applied outside the nodal points and the lower 
within. Let 


r 2 = 4M 

so that eq. (^5) becomes 


cos 0 + sin © - 

r d© 2 
r 


cos © 


sin © - 

r d© 




kl 


where r is now measured in units of r . We have made use of the fact 

o 

that 


, . 1 dr 

tan if — — 

Y r d0 


Between the nodal points the solution is 


r as 1 


and outside the nodal points 


r sin 0=2 


sin 0 
r 


This boundary has been plotted in Fig. 12. 

It is interesting to note that if there were equal plasma winds 
streaming from the right and left the approximate solution would be a 
circle of radius r=l. This is also the exact solution. (The reason 
for the agreement is that the field of a two dimensional dipole confined 
to a cylinder has a magnitude on the surface which is twice the tangential 
component of the dipole alone. This is just Beard f s assumption.) If 
the wind from the left is cut off, the magnetic field will expand into 
this region and the magnetic pressure on the right will drop. The 
cavity surface on the right will therefore fall toward t&e dipole and 
reach equilibrium at the configuration given by the exact boundary 
plotted in Fig. 12. 
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STABILITY 

The stability of the cavity surface is a very important question 

and also a very difficult one. If one is to make any progress at all with 

the analysis one is forced to make some rather drastic assumptions. 

Considerable effort has been devoted to problems of stability. 

Before proceeding with the special features of our problem we shall 

describe an important stability requirement for thermal plasmas 

confined by a magnetic field, (in our problem there is a net mass 

flow In the plasma wind. We are now considering a plasma at rest.) 

This problem has been investigated by Kruskal and Schwarz schild\ 

2 3 

Rosenbluth , Grad and many others. One common feature of these 
treatments is that the surface is stable when the center of curvature 
of the boundary lies in the magnetic field and unstable when it lies 
in the plasma . We may give qualitative evidence for this conclusion 
In this way. We consider first the case where the center of 
curvature lies in the plasma, (see Fig. 21) The magnetic field Is 
directed into the page. In the equilibrium configuration (solid line) 
the magnetic pressure just balances the plasma pressure, i.e. 

B / 8jt = p. The magnetic field Increases as we approach the boundary. 

Thus If the boundary is perturbed as indicated by the dotted line the 
field at b will increase and the field at a decrease. If we adjust 
the deformation so that the plasma volume is kept constant the plasma 

1 M. Druskal and M. Schwarz schild, Proc . Roy. Soc. A223, 3^8, 195^- 

2 M. N. Rosenbluth and C. L. Longmire, Annals of Phys. 1, 120, 1957* 

^ H. Grad, 2d Geneva Conf., 31:190. 
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pressure will remain the same. It is clear then that the perturbation 
will grow. Let us now consider the other case (Fig. 22), The field 
now increases as we move away from the boundary so that in the perturbed 
state the field pressure is increased at a and decreased at b. Thus 
the surface is stable. 

Since it is the latter case that prevails in our problem we might 
expect the surface to be stable. There is however a difference. Our 
plasma is a wind and not stationary. Further the mean free paths are 
too long to apply the continuum plasma equations. Also the magnetic 
field is itself imbedded in a plasma. (The earth* s atmosphere Is 
certainly highly ionized at 9 earth 1 s radii.) 

Parker* 1 " has considered the stability of the interface between the 
solar wind and the earth’s field and concludes that It Is everywhere 
unstable. We do not feel that this conclusion necessarily follows 
from his analysis. We shall conclude now with a more detailed consider- 
ation of Parker’s model. 

The cavity surface is assumed to be locally plane with a plasma 
stream incident on it at an angle 0 . See Fig. 23. This system is 
assumed to be in equilibrium with the magnetic pressure of the wind. 

The magnetic field below the interface is immersed In an infinitely 
conducting. Incompressible, invisid plasma. It Is further assumed that 
this plasma has a scalar pressure and scalar conductivity. The surface 
is now perturbed by a traveling wave. The equation of the surface is 

^ Op. cit. 
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z a r|(y,t) =* Ae^^ 4 * 1 ^ . Tiie equations of motion together with the 
■boundary conditions determine a dispersion relation to = to(k). If a 
solution exists for to complex the surface will be unstable if Ira 
is negative since e ico ^ will now grow exponentially. Parker states 
that "the field density B remains uniform" under such a perturbation. 

This would be true only if the field is directed into the page. Although 
such a relative orientation of field lines and wind direction does occur 
In the geomagnetic problem (in the equitorial plane) it does not occur in 
our two dimensional problem. However, even In the equitorial plane of 
the geomagnetic field, we would suspect that the curvature of the boundary 
would have a considerable stabilizing influence. 

We shall show now that even if one accepts Parker *s model the cavity 
surface surrounding the two dimensional dipole will for the most part be 
stable. We shall show that there is a critical angle of incidence above 
which the surface is stable and below which the surface is unstable. 

Because the magnetic field in our two dimensional problem lies in 
the plane of incidence of the streaming plasma particles we must 
generalize the geometry considered by Parker. In his analysis the field 
is perpendicular to the plane of incidence. We shall investigate the 
following problem, (see Fig. 23 ) The equilibrium boundary lies in the 
x-y plane. The Incident plasma particles have polar angles © q , <( . 

We assume that the surface is perturbed by a traveling wave moving In 
the y-direction. The equation of the perturbed surface Is z = T](y,t). 

We shall first consider the case where the magnetic field is perpendicular 
to the direction of propagation (in the x-direction) and secondly parallel 
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to the direction of propagation (y-direction) . 

Parker has assumed that the plasma beneath the boundary is 

incompressible j invisid, and infinitely conducting. In the unperturbed 

state u = 0, p = P and B = B • In the perturbed state ve assume that 
7 o o 

all of these quantities differ only slightly from their equilibrium values, 
Let the Increments be denoted by n, 5p, and 8B* We shall also find it 
convenient to introduce a displacement vector 5 such that u = . 

The equations that apply in such a plasma are^ 


P ^ = ~^P + (wS)x? 


^ = ^<(ux3) 


V*u = 0 


V-B = 0 


where p is the plasma density, u the fluid velocity, p the plasma pressure 
and d/d t the mobile operator ^ + u*V. Neglecting terms of second order 
the plasma equations become 


P il = _V5 P + jL (Vk5B) X l c 

Bt 


SB = Vx(lx? o ) 


v.g = o 

V. 8B »= 0 


S. Chandrasekhar^ Plasma Physics (University of Chicago Press; i960) . 
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Expanding eq. ( 47 ) we find that 


6B = (i o -v)t - B o (v*t) . 




Since B q has only a x-component 


6B - B 

o ox 


But t =T ( y,z,t ) so that 


6B = 0 

The magnetic field is unchanged to first order. 

It follows then from eq. (^6) by taking the curl of both sides that 


Vx T = 0 

We see then that the fluid motion is irrotational. We may therefore 
represent the fluid motion by a scalar potential such that u « -TO. 

Eq. (^5) becomes 

- pMJ> ss - V&p 
Integrating 

p4> = 6p. (50) 

Also since the motion is incompressible 

V 2 * = 0 . (51) 
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Before one can complete the solution of eqs, (50) and ( 51 ) for the 
plasma motion one must specify the boundary conditions. The boundary 
in our problem is a free boundary and is determined by the condition 
that the pressure on the boundary from the plasma wind must balance 
the pressure of the plasma below. (The pressure from the plasma below 
the boundary is in part due to the particle pressure and in part due 
to the magnetic field pressure.) We must now investigate the effect 
of the perturbation on the plasma wind pressure. 

We consider a small segment of the perturbed boundary as illustrated 
in Fig. 24. The normal to the perturbed surface makes an angle a(y,t) 
with the z-axis. We assume that a is small. Let U be the normal 
component of the velocity of the surface. Parker points out that in 
the local frame of reference moving with the fluid the electric field will 
vanish so that in this frame the particles are specularly reflected. 

The pressure then of this moving surface will be 

p = 2Mm(V n -u) 2 

where is the normal component of the particle velocity and so 

- U is the velocity relative to the surface. The change in pressure 
then is 


6p = p-p Q = 2Nm 

(V -u) 2 -V 2 
v n n 

s 4lfa 

V (V -V )-V u 
n v n n ' n 


__ 0 


_ o 0 0 _ 


to first order in small quantities, V is the normal velocity component 

o 

to the unperturbed surface. The change in the wind pressure is then 
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6p = 4NmV fv cos © sin 0 sin & sin a + U cos 0 1 

* L o o o °A 

and since the magnetic pressure does not change this must equal the 
change in the plasma pressure "below the boundary as given by eq. (50). 
We now let 

Sp = A ± e i(^t+ky) + k Z 

$ = A 2 e i ^ t+ky ^ +kz 

and 

i(r, t) - a 3 e 1 ^^) 

The z-dependence e^ Z has been choosen so that V 2 ® = 0. The three 
constants are determined by two boundary conditions and the equation of 
motion of the fluid. The two boundary conditions are that the pressure 
below the surface must equal the pressure above and that the boundary 
must move with the fluid, i.e. 



on the surface. The equation of motion (eq. (49)) gives 
A 1 /p = i A 2 


53 



k9 


and the houndary conditions become 


itoA^ » -kA 0 and p/L » 4NmV^ -ik sin © cos © sin d A_ 
0 d 1 L O O o l 


k cos © 


V 2 


.] 


( 52 ) 


where we have used the fact that 




The determinant of the coefficients must he zero so that 


p P 

2 4Vk » 4v k 

i ur + cos 0 co 4 

€ O € 


sin © cos © sin d =0 
o o r o 


where we have set the ratio of the densities p/ NM = €• Solving for 
ico ve find 


4 ZVk 

ico = cos 0 

€ O 


1 +- 7 : 

L 


2 x 2^ J 2 j 
€ tan © sin 6 +1 

r\ ~ 


-t- VyT7 




€ 2 tan 2 © sin^d - 1 
o r o 


We see that ico will have a positive real part unless 


e tan 0 sin d = 0 
o *o 
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or 

cos 0=0 
o 

So that the surface wave will grow except when 



or 

o 

Fortunately ^ = 0 or it in both our two dimensional problems , so that 
the surface is stable with respect to this mode. 

Let us now turn to the mode which propagates along the magnetic 

“+ A 

field lines. We therefore take the field to be B = B y, The change 

o 

in the magnetic field 5B is no longer zero. From eq. (49) 



5B = (B o .V)T = B o g (y,z,t) 

Now 


Vx5B = B 

o oy 


0 * 1 ) 


and 


(VxSB)XB q 



-B J- 
o dy 


V(B -1) - B 
' o ' o 
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v X 1 = 0 


2 
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Eq. (55) is Just the dispersion relation for an Alfven wave. If 

such a solution exists, i.e. If we can satisfy the boundary conditions 

with such a mode, it will be stable since w is real. As we are looking 

for possible unstable modes we shall take V x | = 0, so that we are again 

dealing with irrotational flow. The situation is then as before except 
P 2 

that (B“ - B )/ 8 jt, the change in the magnetic pressure, is not zero. 

Here 


= (B -V)| = B 
x o ' o 




so that 


(B 2 -B 2 ) B -5B B 2 Si B 2 dg 
o = o _ _o_ V = _o y 

8n W kit ’by = 4jtiu by 



4jtiw ^2 


4niw 


In the equation of pressure "balance ve must nov add a change in magnetic 
pressure so that eq. ( 51 ) Becomes 


2 2 
B k 

+ Wiw A 2 


i t-NmV <: 


sin 0 cos © 
o o 


sin & — A - 

r o w 2 


kA, 


cos © 


Eliminating the A's and solving for w ve find 


iu) 


Wk 


cos 

£ 


1 + 


■v /2 
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where 


x ^ € tan © sin d 

l-€ o 

The real part of iu> is positive, and hence the surface is unstable 

when 


V2 



> 1 


or when 


2 

tan © I sin <L I > 

o 1 r o 1 r 
ve 

In our problem < j >^ = tt/2 so that instability occurs when ©^ is greater 
than the critical angle ( © Q ) c where 

tan (© ) = 2 / -Sk 

We might say then that for all angles less than (© Q ) c the surface will 
be stable and for all angles greater than (0 ) the surface will be 
stable or unstable depending on the stabilizing influence of the 
curvature of the surface which has been neglected. 

Unfortunately there is no reliable estimates for the plasma density 
as far out as 9 earth’s radii. We would suspect however that the 
density is quite low so that e is small and the critical angle ~ 90°. 
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On the basis of the above calculation we might draw some tentative 
conclusions concerning the stability of the surface surrounding the 
earth's dipole field. In the azimuthal plane passing through the sun 
the surface will be stable over a wide range. In the equitorial plane 
however we feel that the stability question can not be answered by 
the above analysis. The curvature of the surface will exert a 
stabilizing influence which should not be neglected. 

SUMMARY 

We have shown that a plasma Incident on the magnetic field of 
an infinite line current or a two dimensional dipole will confine 
the field so that it forms a cavity In the wind. The plasma particles 
move in straight lines (assuming long mean free paths) and penetrate 
a short distance into the field to form a current sheath within which 
they are deflected out of the field. We have shown that the thickness 
of this sheath Is small compared to the dimensions of the cavity. 

In the steady state the particles will be specularly reflected 
from the surface. They will therefore exert a pressure on the surface 
which is proportional to the square of the normal velocity component. 
The cavity boundary is determined then by the condition that the 
magnetic pressure from within the cavity balance the particle pressure 
from without. We have calculated the shape of this boundary separating 
the plasma wind from the magnetic field of an Infinite line current and 
a two dimensional dipole. The analysis also determines the magnetic 


field within the cavity. 
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The cavity boundary has been shown to be stable except for those 
regions where the cavity wall is almost parallel to the plasma wind 
direction. 
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Fig. 10 Dipole cavity boundary mapped onto 
a circle. The dotted lines are 
branch cuts. 
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Fig. 12 Exact cavity surface as determined from Eq. (4l) and Beard's approximation. 
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Fig. 14 The effect of the field due to the current 
sheath (B) will be to Increase the dipole 
field on the equator and decrease it at 
the poles . 





Pig. 16 Definition of polar vectora used In 

Eq. ( 44 ). The dotted lines are branch 
cuts. varies between - d and + d -it, 

between -d and it - d and^ between 
-tt/2 and -fir/ 2 . 












Fig. 21. Perturbation of equilibrium boundary surface 

when the center of curviture lies in the plasma 
domain. 
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Fig. 22 Perturbation of equilibrium boundary surface 

when the center of curviture lies in the field 
domain . 
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Pig. 24 Perturbed element of boundary surface, n is 
a unit vector normal to the perturbed surface 
and n Q a unit vector normal to the unperturbed 
surface . >7 is the displacement of the perturbed 
surface element from the unperturbed surface. 
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